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Abstract 

The effect of sequence heterogeneity on polynucleotide translocation across a 
pore and on simple models of molecular motors such as helicases, DNA poly- 
merase/exonuclease and RNA polymerase is studied in detail. Pore translocation 
of RNA or DNA is biased due to the different chemical environments on the two 
sides of the membrane, while the molecular motor motion is biased through a cou- 
pling to chemical energy. An externally applied force can oppose these biases. 
For both systems we solve lattice models exactly both with and without disorder. 
The models incorporate explicitly the coupling to the different chemical environ- 
ments for polymer translocation and the coupling to the chemical energy (as well 
as nucleotide pairing energies) for molecular motors. Using the exact solutions and 
general arguments we show that the heterogeneity leads to anomalous dynamics. 
Most notably, over a range of forces around the stall force (or stall tension for 
DNA polymerase/exonuclease systems) the displacement grows sublinearly as i M 
with fx < 1. The range over which this behavior can be observed experimentally 
is estimated for several systems and argued to be detectable for appropriate forces 
and buffers. Similar sequence heterogeneity effects may arise in the packing of viral 
DNA. 
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I. INTRODUCTION 



The dynamics of many single molecule experiments can be described in terms of a 
"particle" moving along a one-dimensional substrate. For example, polymer transloca- 
tion through a narrow pore can be parameterized by the number of monomers threaded 
through the pore. The motion of molecular motors such as kinesins, dyneins, myosin, he- 
licase, DNA polymerase, exonuclease and RNA polymerase can be described by the loca- 
tion of the motor on the one-dimensional substrate (microtubules, actin filaments, DNA 
and mRNA) on which they move. Similarly, the packing of a newly replicated DNA or 
RNA in viruses may be described by the molecular weight of th e pack e d genome. These 
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Under most conditions the motion of the coordinate describing the system is biased in one 
direction. The bias in the case of molecular motors and packing of newly replicated viral 
genomes is due to a chemical process such as ATP (or more generally, NTP) hydrolysis, 
while for polymer translocation it can be generated by the different chemical environments 



on the two sides of the pore. For tran slocating sing! 



provided by adding, for example RecA (jHegner et al 



e stranded DNA, such a bias could be 



19991 ) or other single stranded binding 



proteins (which do not pass through the pore) to the solution on one side of the membrane. 
Single molecule experiments allow another source of bias to be introduced into the system, 
namely an externally applied force F . This has been done, for example, by attaching a bead 



to a molecul ar motor ()V 



the viruses f Smith et al. 



sschcr 



et al 



19991 1 or to the end of the genome which is packed into 

arly, charged 



20011 ) and pulling on it using optical tweezers. Simi 



polym ers have been translocated using an externally applied electric field ((Meller et al 



20011 ). An in teresting var i ant on these experiments is the single molecule measurements of 



Wuite et al. (TWuite et al 



200(1 ) on DNA polymerase, which converts NTPs (nucleotide tri- 
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phosp hates) into a ligated chain of nucleotides via complementary base pairing ((Maier et al 



20001 ) . Wuite et al. apply a force F' not to the motor itself, but instead across the ends of 



the ssDNA/dsDNA complex to create a tension across the substrate on which the molecular 
machine operates. Beyond a critical tension F' c of order of 40 pN, the motor goes backward 
and turns into an exonuclease. The severe stretching of the backbone of the complementary 
DNA strand for F' > F' c presumably makes further conversion of NTPs unfavorable and 
causes removal of nucleotides by the motor to be favored. Forward and reverse motion of 
this enzyme are believed to be associated with different active sites (jDoublie et all 119981 ). 

Most theoretical treatments of these systems have assumed homogeneous (or at least 
periodic) systems. Independent of the microscopic details, such problems can be described 
at long times by a random walker moving along a tilted potential or, equivalently, a biased 
random walker. For molecular motors such as kinesins, dyneins or myosins the assumption 
of homogeneity is indeed, in most experiments, entirely appropriate. However, in other cases 
the motion is along a one- dimensional disordered substrate. This is the case, for example, for 
RNA polymerases, exonuclease and DNA polymerases, helicases, the motion of ribosomes 
along mRNA, the translocation of RNA or DNA through a pore, and the packing of a 
viral genome. In all these systems the one dimensional substrate reflects the heterogeneity 
of DNA or RNA, and leads to a modification of the coarse-grained effective potential in 
which the random walker describing the system moves. The potential now depends in a 
complicated way on the location along the substrate. Two examples of potential energy 
landscapes of particular interest to us here are random energy and random forcing energy 
landscapes. We define a random energy landscape to be any effectively one-dimensional 
potential with a mean slope and fluctuations in the value of the potential with a finite 
variance about this linear tilt. A random forcing energy landscape has an overall mean 
slope but with energy fluctuations which are themselves described by a random walk. In 
this case, the energy fluctuations about a linear tilt grow as the square root of the distance 



along the substrate. These two types of energy landscapes have 



statistical mechanics literature ( Bouchaud et al 



1990; 



jeen studied in detail in the 



Derrida . 



19831 ) and lead to strikingly 



different long time dynamics. In particular the random forcing energy landscape leads to 
behavior quite different from diffusion with drift when the overall tilt of the landscape is 
small, as discussed in detail below. 

Recently, the effect of disorder in the form of defect sites in a rachet model which locally 
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reverse the bias o 



the methods of ( Jiilicher et al 



molecular mo t ors ha s been considered ((Harms and Lipowskvl . 119971 ) , using 



1997). It was suggested that even though fluctuations in 



the microscopic potential are bounded, the resulting effective energy landscape is random 
forcing. Specifically, it was argued that when the defect concentra tion was large enough 
anom alous random force dynamics would appear. As pointed out in (jLubenskv and Nelson , 



20021 ) heterogeneity in base pairing energies also leads to a random force landscape in the 



context of DNA unzipping. 

In this paper we study the effect of sequence heterogeneity in both polymer translocation 
and molecular motors in detail for an exactly solvable class of simple lattice models. We 
consider both systems in the context of single molecule experiments that apply an external 
force pulling back on the polymer or the motor, which in the absence of this force are 
biased to move in one direction. We introduce microscopic models for both systems which 
can be solved exactly both with and without disorder. A generalization of our motor model, 
discussed in Section V and Appendix D, can als o be used as a very simple model of the DNA 



polymerase/exonuclease experiments of Ref. (jWuite et all 120001 ) . One can also consider 
closely related models of the packing of a viral genome. In this case there is an extra source 
of bias due to the energetic cost of packing the DNA inside the virus. The externally applied 
force acts in conjunction with this bias while the motor acts against both. The details are 
very similar to the cases discussed here, with the exception that the energy cost of forcing 
the DNA into the capsid does not necessarily vary strictly linearly with the amount of DNA 
that has entered. We do not include a separate discussion of this interesting system. 

We show that sequence heterogeneity of single stranded DNA or RNA and heterogeneous 
base pairing energies have a dramatic effect on the dynamics of both systems. For a homo- 
geneous substrate and no chemical bias, the average velocity changes monotonically through 
zero as the external force is varied, changing sign as the force reverses direction (see Fig. 
ma)). When a chemical bias (which we take to act in the direction opposing the force) is 
present, the scenario is similar with the velocity changing sign at a stall force, F s , which 
depends on the degree of chemical bias (see Fig. IHb)). In contrast, the combination of a 
disordered substrate and a chemical bias produces very different behavior for both systems. 
In this case we show that generically disorder introduces a random forcing effective energy 



landscape which is res ponsible for the anoma 
Harms and Lipowsky ((Harms and Lipowskvl . 



ous d ynamics. Similar to the observation of 



19971 ) a random forcing landscape is gener- 
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FIG. 1 Schematic behavior of the drift velocity at long times for homogeneous and heterogeneous 
systems as a function of the applied force, where a positive force resists the chemically favored 
direction of motion. It is assumed that chemical forces (such as ATP hydrolysis or chemical binding 
on one side of a pore) lead to a positive velocity in the absence of a force, (a) No externally applied 
chemical bias (Afi = 0). (b) A finite chemical bias (A/i > 0), where the light line corresponds 
to homogeneous or periodic environments and the solid line refers to heterogeneous environments. 
The anomalous dynamics ((x(t)) ~ t^, with < 1) arises in the vicinity of what would be the stall 
force, F s , for the homogeneous system. For F< < F < F>, the effective velocity depends on the 
width of the time averaging window, and tends to zero as the width of the window goes to infinity. 
The striped line denotes the region where anomalous diffusion is also present. 



ated even if we neglect an explicit contribution ((Lubenskv and Nelsonl . I2002T ) from random 
base pairing energies. We discuss three different dynamical regimes which arise due to this 
landscape as the externally applied force is varied. The most notable transition arises in 
the velocity of the random walker describing the system. Specifically, we find that there 
are critical values of the force F> and F< such that for any force between these values the 
velocity is zero in the sense that the average particle position (x(t)), where (. . .) denotes an 
average over thermal fluctuations, increases as a sublinear power of time. We also discuss an 
even broader range of forces where the diffusion is anomalous (see Fig. IHb)). The transition 
points between the different types of long-time dynamics can be calculated exactly for the 
simple models studied here. 

Under special conditions a random energy landscape is also possible. In this case the 
expected behavior as a function of force is similar to a homogeneous system: The potential 
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FIG. 2 Schematic behavior of velocity for homogeneous and heterogeneous systems as a function 
of the chemical bias Afi. (a) No externally applied force, (b) A finite externally applied force 
for homogeneous (light line) and heterogeneous (solid line) substrates. The anomalous dynamics 
arises in the vicinity of what would be the stall chemical bias, A/x s , for the homogeneous system. 
As in Fig. ^ the striped line denotes the region where anomalous diffusion is present. 

fluctuations simply renormalize the drift velocity and diffusion constant at long times. That 
is, as the applied force is varied the behavior is similar to that of a homogeneous system with 
no chemical bias. Provided that random contributions to the energy landscape not associated 
with simple conversion of chemical energy can be neglected, random energy models describe 
the dynamics in the absence of chemical bias (see Fig. CJa)) on heterogeneous substrates. 

An alternative way to observe the anomalous dynamics is by holding the external force 
constant and varying the chemical bias. This can be done by changing the concentration 
of, say, nucleotide triphosphates for molecular motors, or by changing the concentration of 
the polymer binding protein in one chamber for polymer translocation experiments. In this 
case, when the force is held at zero, the velocity changes monotonically in tandem with 
the chemical bias (see Fig. I2fa)). However, when the external force is held constant at a 
non-zero value, a region with anomalous dynamics appears as the chemical bias is varied 
(see Fig. Efb)). Between two values of the chemical bias A/i< and A/i>, the displacement 
of the particle with time is again sublinear, in contrast to the same experiment performed 
on a homogeneous substrate. As illustrated in Fig. [21(b) the velocity is then a monotonic 
function of the chemical bias, changing sign at a stalling chemical bias Afi s . A summary of 
the qualitative behavior of the velocity as function of both the chemical bias Afi and the 
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FIG. 3 The dependence of the velocity on the chemical bias A/x and the external force F. We 
neglect for simplicity contributions to a random force landscape (such as fluctuations in base pairing 
energies) which may be present even for A// = 0. Here it is assumed that the chemical bias always 
act in the direction opposing the force. The black wedge denotes a region of sublinear drift with 
time, i.e., effectively zero velocity. 

external force F is shown in Fig. El It is worth noting that there is no region of sublinear 
displacement when A/i = because the energy landscape is then random energy rather than 
random forcing, whereas when F — 0, there is still a random forcing landscape everywhere 
except exactly at stalling, but the randomness is too small in the vicinity of A/i = to cause 
anomalous dynamics. 

To keep the discussion simple, Fig. Elneglects contributions to a random forcing landscape 
other than those produced by the simple conversion of chemical energy along an inhomo- 
geneous track. Additional random forcing contributions will arise from, e.g., base pairing 
energies in the case of helicases, which open up DNA strands or DNA polymerases and 
exonucleases, which add or delete complementary base pairs. Motors, such as RNA poly- 
merases and ribosomes produce trailing strands of mRNA and protein respectively. Since 
these products are themselves heteropolymers, composed of monomers which interact dif- 
ferently with the solvent, here too we would expect additional contributions to a random 
forcing landscape. Such effects will only accentuate the anomalous dynamics which is the 
subject of this paper. 

Before concluding this introduction, we should emphasize our perspective on the models 
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of polynucleotide translocation and molecular motors studied here. In an effort to obtain 
simple, soluble models which incorporate heterogeneity, we intentionally neglect important 
molecular details such as those which describe the detailed pore interactions of the translo- 
cating nucleotides or distinguish the biological role of motors such as helicases, DNA poly- 
merase and exonucleases, RNA polymerase, etc. The motors mentioned above perform 
important specialized functions such as opening double stranded DNA, polymerization and 
depolymerization or creating messenger RNA while moving along heterogeneous tracks. Such 
functions are incorporated into our model simply by adding an explicit (position-dependent) 
chemical force to the e nergy landscape. More sophistica t ed at tempt s to get molecular details 
right (see, e.g., Refs. (jGoel et all 120031 ) . (pinion et all Il992r ) and ([Betterton and Jiilicher . 



20031 )) serve a valuable purpose, which can be important for modelling some aspects of 



the dynamics on various time scales. However, incorporation of sequence heterogeneity, ne- 
glected in most previous modelling efforts, is nevertheless crucial to correctly describe the 
anomalous long time dynamics (e.g., (x(t)) ~ t M with \i < 1) near the stall forces in these 
systems. Otherwise, we expect simple diffusion with drift (similar to what we find here for 
homogeneous models or a random energy landscape) at long times. We do not expect the 
multiple intermediate states and numerous rate constants of more sophisticated models to 
change our predictions of heterogeneity-induced anomalous dynamics at long times. 

The paper is organized as follows: In Section II, to establish notation and provide a con- 
text for the rest of the paper, we discuss the homogeneous models for polymer translocation 
and molecular motors is some detail. In Section III the effect of the heterogeneity on the 
energy landscape is introduced. Section IV discusses the resulting dynamical behavior and 
the exact location of the transition points within the models. Finally, Section V estimates 
the experimental range over which the anomalous dynamics may be observed for a few repre- 
sentative biological systems and discusses the effect of finite time experiments on the shape 
of the velocity-force curve. 



II. HOMOGENEOUS MODELS 



Before turning to heterogeneous systems we first define microscopic models for both 
homogeneous polymer translocation and molecular motors. The simplicity of both models 
allows for their exact solution. Dynamics in heterogenous systems will be treated in Sections 
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FIG. 4 Schematic picture of the polymer translocation experimental setup considered. A polymer 
is biased to move through the pore by a solution of binding proteins in the right chamber. A bead 
exerts a force in the opposite direction. The arrows reflect the lack of inversion symmetry in, e.g., 
single stranded DNA or RNA. 

Ill and IV. 



A. Polymer Translocation 



An idealized experimental setup is shown schematically in Fig. EJ A polymer is threading 
through a narrow pore located on a two dimensional membrane which separates two chemi- 
cally distinct solutions. For concreteness we consider the right side as containing a polymer 
binding protein which is absent in the left hand side. In addition, a bead, through which 
a resisting force is exerted on the polymer, is connec ted to the lef t end of the polymer. A 
model of this kind has been discussed by P. Nelson (jNelsonl . 120031) as a simple ex ample of 



stochastic rachet-like dynamics in biological systems (see also (jPeskin et all 1199 3)). Alter 



natively a force c ould be applied via an ex ternal electric field acting across the pore on a 



charged polymer (jKasianowicz et al 



19961 ) 



A convenient representation of the system is through a one-dimensional random walker 
located at a coordinate x which represents the length of the polymer that has translocated 
to the righthand side. The conditions under which the full three-dimensional, multi-species 
problem can be simplified are reviewed below. The dynamics of the random walker is 
governed by the interaction of the polymer with the pore, the binding of the protein in the 
right chamber, and the externally applied force. 

Before turning to a specific microscopic model consider the general form of the potential 
experienced by the random walker due to all these interactions. Because we neglect sequence 
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FIG. 5 (a) The periodic potential due to pore interactions with a translocating polymer without 
inversion symmetry, (b) The tilt of this potential generated by a combination with the binding 
protein and the external force. 

heterogeneity in this section, the energy due to interactions with the pore, U(x), is some 
periodic function with a period given by the size of a monomer. An example is the sawtooth 
or rachet potential shown in^a). This type of potential accounts for an energetic barrier 
for translocation through the pore. The lack of inversion symmetry reflects, for example, 
the difference in passing single stranded DNA or RNA in the 3' — ► 5' direction through 
the pore as opposed to the reverse. The energy due to the interaction with the polymer 
binding protein is however very different and has the form —F^x, growing linearly with x. 
Thus the energy decreases as the polymer translocates to the righthand side. The value of 
is governed by the chemical potential difference per monomer, A/i, of the polymer in 
the solutions on the righthand and lefthand sides. This chemical potential difference is a 
function of the protein concentration and its binding energy to the polymer (a more detailed 
description of F^ for the microscopic model discussed below is presented in Appendix IA"|) . 
Finally, the backward force applied on the bead leads to a contribution to the energy of the 
form Fx. Upon collecting together these contributions, the total potential experienced by 
the random walker, is given by 



As is evident from the effective energy landscape shown in Fig. E[b), the random walker is 
moving in a periodic potential with an overall slope which depends on the protein concen- 
tration and binding energy as well as the external force. Such a potential leads on long time 
scale and large length scales to motion which is diffusion superimposed on an overall drift 
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FIG. 6 Graphical representation of a simplified model for polymer translocation or molecular 
motors. These two cases are distinguished by the choice of rate constants (see text). The distinct 
even and odd sublattices are denoted by a and b respectively. 

velocity. Thus, the average location of the particle (x) behaves as (x) = vt while the mean 
square fluctuations about this drift behave as (x 2 ) — (x) 2 = 2Dt, where v and D depen d 
on — F and the details of the rachet potential (see, e.g., (jLubenskv and Nelsonl . Il999j )). 
Here, the brackets (. . .) represent an average over thermal fluctuations. 

We emphasize that our simplified description in terms of a single coordinate x that 
diffuses and drifts in a one-dimensional energy landscape is valid only when the transla- 
tion al motion of the polymer bac kbone through the pore is the slowest process in the prob- 



lem rtLubenskv and Nelson 



1999). In particular, this model assumes that the translocating 



polymer is not so long that the relaxation times in the cis (left) or trans (right) chambers 
exceed the diffusion time for the backbone through the pore. This simplified model is also 
inadequate if the polymer can become bound to the pore interior for long perio ds, as recent 



exper iments suggest occurs for one of the best studied polymer-pore systems (jBates et al 



20031 ) . In this case, x will still undergo biased diffusion on long enough time scales, but its 
velocity and diffusion coefficient will no longer be determined by a simple potential U(x). 
Finally, the effect of binding proteins can be captured by a single free energy parameter 
A/x only when their binding and unbinding kinetics are sufficiently fast. The opposite limit, 
in w hich proteins bind irreversibly, but slowly to the polymer, has also received atten- 



tion rtPeskin et al 



1993; 



Simon et al 



199 



2: 



Sung and Park. 



1996I ). but we will not consider 



it further here. 

We now define a simplified microscopic model for the motion of a random walker 



in such a potential. Our mode 



in ( Fisher and Kolomeiskv 



1999; 



is in the spirit of thos e ana 



Kolomeiskv and Fisher 



yzed for motor proteins 



19991 ) (see also Sec. II.B), 
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and allows exact result s for t he diffusion and drift on lon g times. In the language of 



( Fisher and Kolomeiskv . 



1999; 



Kolomeiskv and Fisher . 



19991 ) our model is a n = 2 model 



corresponding to a motor with just two internal states. More importantly, our model gener- 
alizes naturally to a heterogeneous version (see Section III) for which exact results are also 
possible. We allow x to assume a discrete set of values x m , where m = 0,1,2... labels 
distinct a (even) and b (odd) sites. We can allow different distances between x m+ \ — x m 
and x m+ 2 — x m+ i but require x m+ 2 — x m = 2a>o which we assume for simplicity is the size of 
the polymer unit which accommodates a single adsorbed protein. For a homopolymer the 
interactions with the pore are some periodic function with a period which we can take to be 
2ao- To model this situation we take odd labelled sites to have a higher energy than even 
labelled sites. The arrangement is shown schematically in Fig. El Even sites have an energy 
e = while odd sites (corresponding roughly to the peaks in the rachet potential of Fig. |SJ) 
have a higher energy e = Ae. Also, indicated in the figure are the hopping rates which de- 
scribe the dynamics of the random walker. The detailed balance condition (in temperature 
units such that ks = 1) is satisfied by 



K 
K 



toe 



-Ae/T-f/2T 



J/2T 



c y e (-A A1 -A £ +//2)/T 



(2) 



uj e 



-f/zr 



Because of the lack of reflection symmetry in the translocating DNA or RNA (for our model 
this asymmetry could be represented by taking x\ — Xq ^ X2 — £i), we expect the intrinsic 
hopping rates to be unequal, uj ^ uj' . The bias induced by the interaction of individual 
monomers with the reservoir of proteins on one side of the pore has been accounted for by 
the chemical potential difference A/i. A more detailed discussion of the dependence of A/x 
on the protein binding energy and its concentration is given in Appendix EI The effect of the 
applied force is included through the parameter / = Fclq. Note that the bias controlled by 
A/x > arises only for steps from odd t o even sites since a prote i n is a ssumed to bind only to 
a whole monomer. As pointed out, in (jKolomeiskv and Fisherl . 119991 1 . other /-dependences 
of the rates consistent with detailed balance are possible. We shall be content with the 
simple one displayed in Eq. El which corresponds to choosing x± — Xq = x^ — x\. 

To show that this microscopic model embodies an effective potential of the form (JTJ, we 
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eliminate the odd-numbered sites. This elimination can be accomplished by formally solving 
the equations of motion for the odd sites, substituting into the remaining even site equations, 
and taking he long-time limit (see Appendix^. Alternatively we can invoke detailed balance 
and consider an effective energy difference AE = E(m + 2) — E(m) between site m + 2 and 
m where m is even. Upon setting 



W, 



m,m+2 _ e (E(m+2)-E{m))/T 



W m+2 

where W„ m is the effective transition rate between site m and n, we have 



(3) 



AE 



E(m + 2) - E(m) 



(4) 



Use of the rates Q leads to 

AE = -A/2 + 2f (5) 

as one would expect. Note that when the force vanishes (/ = 0) and the chemical potential 
gradient A/i = one has AE = and no net motion is generated. More generally, an 
effective tilted potential of the form is generated, with A/x > causing a drift of the 
polymer to the right. The external force on the left can reduce or even reverse the overall 
slope. Such a potential inserted into microscopic rate equations for the even sites (see 
Appendix |Bj) is well known to lead to diffusion with drift on long time scales and large 



In fact, for this 


mode 


using the results of ( 


Derrida 


Fisher and Kolomeiskv 




1999; 


Kolomeiskv and Fisher 




1999 


) one 



one can calculate the veloc- 



ity and diffusion constant exactly. After some lengthy calculations, one obtains for the 
velocity 

v = 2a - 



w„ wi 



w^ + w^ + w b +w b 



The diffusion constant of the model is given by 



with 



D = 2a 2 (wT^r + w ^ w T) + ^w^w^w^w^ K 

Or + w t + w 7 + w v) 3 



K = (w-f + (wD 2 + «) 2 + «) 2 

+ 2 i. w 7 w 7 + w7 w 7 + W 7 W 7 + W 7K 



(6) 



(7) 



(8) 
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It is interesting to set / = and consider the limit of A/x/T <C 1 (small chemical bias, no 
external force) and the limit A/x/T — > oo and A/x ^> Ae (large chemical bias, no external 
force). When A/x/T <C 1 the velocity takes the linear response form 

2a aWe- Ae/T A/x 



In the limit of A/x/T — > oo and A/x ^> Ae, the velocity saturates at f max , with 



(9) 



COCO 

"max = ^ : Af/TV i A ' V"*-^/ 

In both cases the velocity is a decreasing function of Ae, as one might expect because the 
rate limiting step in this simple polymer translocation model is the energetic barrier as each 
successive segment passes through the pore potential. 

For the diffusion constant one finds similarly in the limit A/x/T <C 1 

2 uju'e- A£ / T 

2 V MHe-^)+,'(l-e- A£/T )) nn 

Like the velocity, in this regime the diffusion constant decreases as Ae increases. Note that 
the diffusion constant deceases when A/x increases. This behavior arises since the rate of 
backward steps decreases as A/x increases. In the limit A/x/T — > oo, and A/x ^> Ae we find 
that the diffusion constant saturates at 

2 W«a; + oo') 2 e 2A ^ T + u\l + 2e Ag / T )) 
max_a ° {u + e^l T {uj + uj')f ' ( j 

which also decreases with Ae. The diffusion constant again decreases as a function of Ae 
due to the rate limiting step of the passage through the pore. 



B. Molecular Motors 



A typical experimental setup is shown in Fig. The motor attempts to move from 
the + end to the — end by utilizing the chemical energy stored in ATP or some other 
source of chemical energy. For RNA polymerase, this energy source would be the nucleotide 
triphosphates which are converted into mRNA (not shown). A force (say from an optical 
tweezer) pulls in the opposite direction to the motion generated by the ATP. In this section, 
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FIG. 7 Setup modelled. The motor is moving from the + end to the — end. A force is pulling on 
the motor in the opposite direction. Note that some of the specific biological examples considered 
in the text are more complicated and may be driven by energy sources other than ATP. 

we focus primarily on models of relatively simple motors as in Fig. and mention only in 
passing more complicated effects associated with motors su c h as h elicases or RNAp. 



Theoretical models of molecular motors ( .liilicher et al 



J7) 



1997J) have demonstrated how 



an effective potential of the form (JJ is generated as a result of the coupling to an energy 
source like ATP for a general class of periodic substrate potentials which lack inversion 
symmetry. Here we again introduce a simple model for a two level rachet which is amendable 
to an exact solution, similar t o an n = 2 version of the mod els of Fisher and Kolomeisky 



( Fisher and Kolomeiskv 



199 



Kolomeiskv and Fisher 



19991 ). Like the model for polymer 



translocation in Section A, this motor model will allow us to study the effect of heterogeneity. 
We first consider the homogeneous motor model in some detail. 

We again consider a one-dimensional lattice where even sites have energy e = while 
odd sites have an energy e = Ae. The odd sites represent an "inchworm"-like walking which 
is facilitated by chemical energy released by, e.g., hydrolysis of ATP. The transition rates 
depicted in Fig. El now take a different form, namely 



K 
K 



f/2T 



[a + u))e 



(a'e A ^ T + u')e- A ^ T+ ^ 2T 



(13) 



Note that there are two parallel channels for the transitions (I.Tulicher et all Il997f ). The 
first, represented by contributions containing a and a', arise from utilization of chemical 
energy. The second channel, represented by the terms containing u and u', correspond to 
thermal transitions unassisted by the chemical energy. A/i is given by the standard relation 
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( Howard 1200 lh . 



A/i 



/ [ATP] \ ( [ATP] t 
n \[ADP][P]J n \[ADP] cq [P] 



CC| 



(14) 



where the square brackets [. . .] denote concentrations under experimental conditions and the 
[. . .] denote the corresponding concentrations at equilibrium. We have again assumed the 
external applied force / biases the motion in a particularly simple way. If the substrate lacks 
inversion symmetry, we have a' 7^ a and u/ ^ u. As discussed in the introduction, in some 
cases an additional force arises from, e.g., base pairing energies in the case of helicases, DNA 
polymerase and exonucleases. Similarly an addition force arises also for motors such as RNA 
polymerase and ribosomes which produce trailing strands of mRNA or protein respectively. 
Here we ignore such contributions although they could easily be added in a simple way to 
the model through a redefinition of / through / — > f + where f ^ is the additional force. 
The model is formally similar to the model of polymer translocation, although the different 
functional form of w^, w^~, and w^* has important consequences. 

First we consider the effective energy landscape. To this end, we again eliminate the odd 
sites and describe the remaining dynamics in terms of an effective potential. This is the 
effective potential under which a random walker satisfying detailed balance would exhibit 
the same dynamics. From a formula similar to Eq. 01 ones finds that AE = E(m+2) — E(m), 
where m is an even site, is given by 



AE = Tin 



(15) 



where we have used the rates (fHIj) 

Note that when the external force / = and the ATP/ADP+P chemical potential 
difference A/i = 0, one has AE = and no net motion is generated. Also, when there is 
directional symmetry in the transition rates a = a', uj = oj' and / = one has AE = 
even when A/x 7^ 0. Absent this symmetry, chemical energy can be converted to motion and 
an effective t ilted potential i s generated. Although, t hese conditions are equivalent to those 



presented in (jjiilicher et al. 



1997 



Prost et al. 



19941 ) for continuum models, it is interesting 



to see them at work in the "minima l " mod el studied here (see also (jFisher and Kolomeiskv , 
19991 ) and ((Kolomeiskv and Fisherl . [1999)). The effect of the externally applied force is 
simply to change the overall tilt in the potential. 



1(3 



For a motor on a homogeneous or periodic substrate the effective potential generated by 
the coupling to the chemical potential is thus qualitatively the same as that of a polymer 
translocating through a pore. Again on long time scales and large length scales the dynamics 
is just diffusion with drift. The equation for the velocity and diffusion constant are given 
by Eqs. I6I7I8I together with the rates displayed in Eq. EH 

As for the polymer translocation problem, it is interesting to consider various limits for 
the case / = 0. Using (JT3jl we find in the limit of A/i/T < 1 a drift velocity 

2a (uj'a - uja')e~ Ae/T /Afi\ 



- (16) 



(a + uj + a' + u/)(l + e~ A£ / T ) \ T J 

Therefore, for small A/i/T, the velocity decreases as Ae increases. Note that even when 
A/i ^ 0, v vanishes for a symmetric substrate, i.e. for uj' = uj and a' = a. A natural measure 
of the asymmetry of the potential is u'a/ua'. When this quantity is greater than one (less 
than one) a positive A/i induces a motion to the right (left). This result remains valid to 
any order in A/i. 

The maximum possible motor velocity v max is obtained in the limit A/i/T — > oo and 
A/i ^> Ae, where 



uj'a — Lua' 



w max = 2a ■ — — . (17) 

a + a' 

In contrast to the previous regime and the polymer translocation problem, the velocity is 
insensitive to Ae. Because of the injection of large amounts of external chemical energy, the 
barrier Ae no longer controls a rate limiting step. 

For the diffusion constant of this model of molecular motors in the limit A/i/T C 1 we 
find 

2 (a + u,)(a> + «/)e-**' T 
{a + u + a' W)(l + e- A£ / T ) 

2 A^i 2e~ A£ / T G 

+ a ° T {a + u + a' + cu') 2 (l + e~ A£ / T ) 2 ( } 



with 



G = e- Ae/T (a + u-a[- u')(a'u - auj') 
+ aa'(2(a + a') + 3(uj + uj')) 

+ (u + u')(a'uj + auj') + a' 2 u + ctV (19) 
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Like the velocity, the diffusion constant decreases as Ae increases in this regime. Note that 
the diffusion constant increases as A/i increases, because A/i enhances the rates of motion 
in both directions. In the limit A/i/T — > oo one obtains 



9 ooa' + uo'a + 2aa' 



Anax = 2a 2 — . (20) 

a + a' 

Again, for large chemical potential differences the result is independent of Ae. 



III. THE EFFECT OF HETEROGENEITY ON THE ENERGY LANDSCAPE 

Next we discuss the effect of heterogeneity on the effective energy landscape experienced 
by motors or translocating polymers. The detailed dynamics which results will be consid- 
ered in Sec. IV. As we shall see, heterogeneity has dramatic consequences over a range of 
parameters close to the stall force. 

We first consider the somewhat simpler problem of heterogeneity and polymer translo- 
cation. We then show that a similar picture arises for motor proteins on heterogeneous 
substrates like DNA or RNA. 



A. Polymer Translocation 

Two sources of heterogeneity affect polymer translocation. Both arise for polymers com- 
posed of different types of monomer. We assume for simplicity that the monomers composing 
the polymer are drawn from some random distribution with a finite variance. Provided the 
correlations along the backbone are short range our results are insensitive to the exact na- 
ture of the distribution. The effect of sequence heterogeneity corresponding to a particular 
nucleotide sequence could easily be incorporated into a numerical analysis of the dynamics. 

We first consider general features of the potential for a model with sequence heterogeneity. 
Randomness in the composition of the polymer will of course modify the interaction potential 
between the polymer and pore, U(x). It is easy to see that this leads to a random potential 
component with a finite variance around its mean value, i.e., a random energy landscape. 
The second, more striking effect, arises from the randomness in the binding energy of the 
proteins. The associated force depends specifically on the location x along the polymer. In 
a convenient continuum notation, the total energy gained by attaching to the monomers 



18 



UeffMt 




FIG. 8 Graphical representation of the energy landscape in the case of heterogeneous polymer 
translocation when the chemical environments on both sides of the pore are different. Potential 
fluctuations about the mean slope scale like y/x for large x. The same picture holds for molecular 
motors moving on a heterogeneous substrate powered by a finite chemical potential difference. 

has the form J X Fn(x')dx' , where F^x) represents the different binding energies associated 
with the sequence of the polymer. If the sequence is random the fluctuations around the 
mean slope of the potential grow like \/x. The effective potential experienced by the random 
walker is therefore 

U eS (x) = U(x) - ( f F^x')dx' - Fx) , (21) 

J 

where we have included the externally applied force, F. A schematic representation of 
the potential is shown in Fig. |H1 Since Jq F^x^dx' has fluctuations which grow as yfx, the 
sequential binding of proteins to a translocating polymer creates a random forcing landscape, 
in contrast to the landscape denned by Eq. ^ Because the energy landscape itself can be 
viewed as a simple random walk about a linear landscape, the random force contribution 
to U c fj(x) (an integrated random walk) dominates the random energy term arising from 
interactions with the pore. As will be discussed in Section IV this results in unusual behavior 
if the externally applied force lies in a certain range of values near the stall force. 

Note that it is also possible to obtain a purely random energy landscape in polymer 
translocation. When the chemical environments on both sides match (e.g., for identical 
concentrations of binding proteins) one has F^x) = 0. The only random component of 
the energy landscape is due to the potential for translocating through the pore which has 
bounded fluctuations about its mean value. For this energy landscape, the dynamics at 
long times and large length scales is then biased diffusion, with a drift velo city and diffusion 



constant renormalized by the heterogeneous interactions with the pore (jAlexander et al 
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19811 ). 

We now explore these effects within our microscopic model of polymer translocation. The 
heterogeneity is introduced into the model through the rates (J2J). Imagine drawing the set 
of parameters {p} = {u, u', Ae, Afi} from random distributions (corresponding to various 
nucleotide sequences) with a finite variance. According to Eq. (JHJ),, the total change in 
energy after m monomers translocate is given by 

m 

E(m) = 2/m + Yl AE ( l ) ■ ( 22 ) 
1=1 

Here the AE(m) are effective energy differences between two even sites corresponding to 
the set of values of the set {p} drawn randomly. Since the energy is a sum of independent 
random variables a random forcing landscape is developed. 

We expect that a simple random energy landscape results if we turn off the protein 
binding by setting A/x = 0. However, because the energy at even sites is always E = in 
our simple model, the landscape is just a uniform tilt in this limit. A more realistic model 
would allow additional energy variations at these sites. If we assign an energy e(m) to these 
even sites, it is straightforward to show that the total change in energy after m monomers 
have translocated takes the form 

E(m) = 2/m + e(m) , (23) 
corresponding to a random energy landscape. 

B. Molecular Motors 

We now turn to the effect of heterogeneity on molecular motors. Here, as for polymer 
translocation, we select the set of parameters {p} = {a, a', u, u>', Ae} from a random distri- 
bution with a finite variance. For some motors and enzymes (for example RNA Polymerase, 
helicases, DNA polymerases and exonucleases - see introduction and below), A/i may also be 
random. This clearly only adds an additional contribution to the random forcing landscape. 
Using the results presented above it is easy to see from Eq. (j!5|) that the total effective 
energy change after m monomers is given by 

m 

E{m) = 2/m + ^ AE(l) . (24) 
i=i 
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Here, each AE(m) corresponds to an independent set of values of {p} drawn randomly. 
Thus, as in the polymer translocation problem, the potential is random forcing. 

For motors such as helicases, DNA polymerase and exonucleases, and RNA polymerase 
and ribosomes an additional contribution to the random energy arises due to the force 
associated with, e.g., base pairing energies or the trailing strand which is produced. The 
effect of this would be to modify the expression above to 

m m 

E{m) = 2/m + ]T + E AE ( l ) > ( 25 ) 
1=1 1=1 

where f ^ is the additional contribution of the explicit random forcing from monomer m. 
The resulting random forcing landscape is even more pronounced. 

The above scenario applies as long as the chemical potential difference A/i 7^ 0. In the 
case when A/i = it is easy to see that AE(m) = unless we allow, as in the polymer 
translocation problem, for the energy at even sites also to vary and take the value e(m). In 
this case we obtain 

E{m) = 2fm + e{m) , (26) 

corresponding to a random energy landscape provided e(m) has only short range correlations. 
Although we could write the energy in the form of Eq. now AE(m) = e(m) — e(m — 1), 
so AE(m) is effectively the gradient of a random potential with bounded fluctuations. Note, 
however, that for motors with an f ^ contribution (as in Eq. it is not possible to obtain 

a random energy landscape. 

The energy landscape for both polymer translocation and molecular motors is therefore 
qualitatively identical. Generically, in both random forcing energy landscape de- 

velops. However, if the motor model without the applied external force has no bias (i.e., 
if A/i = 0), we recover the diffusion with drift dynamics associated with a random energy 
potential. 

IV. DYNAMICS IN HETEROGENEOUS ENVIRONMENTS 

In this section we discuss in detail the dynamics of translocating polymers and motor 
proteins with heterogeneity for the model depicted schematically in Fig. |H1 We describe 
four distinct cases with different dynamical behaviors as the externally applied force is 
varied. The critical forces for the transition between the regimes can be calculated exactly 
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in terms of the rates w^*, w£~, w^, averaged over their heterogeneous generalization with 
/ = 0. The explicit expressions for polymer translocation and molecular motors can be 
easily obtained by using the rates in Eqs. (J2J) and (fT3)) respectively We assume throughout 
that A/i ^ 0, as the case A/i = leads only to a random energy model and biased diffusion. 
Also, contributions to the random forcing energy landscape of the form of Eq. (J25|) are 
omitted for simplicity Their addition is straightforward and can be easily seen to enhance 
the region of anomalous dynamics. 

The dynamical behaviors of random walkers in random forcing or ran dom energy land 



scapes have b e en stu died in detail in the statistical mechanics literature (jBouchaud et al 



1990; 



Derridal 119831 ) . Unusual dynamical behavior ar ises for random walkers in a random 



forcing energy landscape. Using the results of Derrida (jDerridal . 119831 ). one can calculate the 
transition points between the different regimes including the effect of randomness. Parts of 
the calculation are outlined in Appendix along with the different regimes in terms of the 
transition rates , , , . Here we consider the experimental setup in Figs. Eland 
E] where the external force is varied. Denoting spatial averages by an overline and using the 
results of the Appendix C one finds the following regimes. 



Regime I: The velocity v and diffusion constant D of the model are finite when 



T 

f < In 

J 4 



wywi 



or 



T 

/> I ln 



(27) 



(28) 



where the subscript / = denotes that / has been set to zero in the average. In this 
regime (x) = vt and (x 2 ) — (x) 2 = 2Dt for long times, where the angular brackets denote an 
average over different thermal histories of the system. Simpler conditions can be obtained 
by assuming that AE(m) = Tin {{w^w^~) / '(w^*w^*)) has a Gaussian distribution about the 

2 

mean 2f + AEj =0 (see Eqs. IS1andlT3j) and a variance V = (AE) 2 j =0 — (AE) f=Q . Here again 
the subscript / = denotes that averages are taken with the value of the force set to zero. 
In the case one has 



/ > \ (AE f=0 + V/T 
V/T 



2 

/ < \ (M /=0 



(29) 
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Note that the force does not contribute to the variance so that V = AE 2 f=o — AE 



7=o 



AEj - AEf . 

Regime II: The velocity v is finite but the diffusion constant is infinite. Thus, in this region 
(x) = vt and (x 2 ) — (x) 2 ~ t 2 ^, where 1 < jjt, < 2. The relevant force ranges are 



-In' 



™7 W b ) /=Q 



■?lnf^ 



™a W b ) f =0 



and 



T ln w^b 



</<^ln 

4 V wTwt 



(30) 



(31) 



J a w b / f=o 

Provided that AE has a Gaussian distribution the conditions reduce to 



\ (AE f=0 + V/2T) < f<\ (AE f=0 + V/T 

V/2T 



-(AE 



/=o 



2 

V^/T) < / < \ (AE f=0 



(32) 



For a Gaussian distribution it is known (jBouchaud et al 

by 



19901 ) that the exponent fi is given 



(33) 



H = 2T\EE f=0 -2f\/V. 

Regime III: The velocity v is zero in the sense that (x) ~ f 1 , where /i < 1. The exponent 
/i also controls the variance, (x 2 ) — (x) 2 ~ t 2tl . This behavior occurs when 



<f<^ln 

™7 w v ) f =o 2 \ w i w rj f =o 



(34) 



When AE has a Gaussian distribution, these conditions reduce to 

- (M /= „ - V/2T) <f<\ (AE f=0 + V/2T) 



(35) 



2 v ' / ' 2 

Sinai diffusion: Here (x) = and (x 2 ) ~ (ln(i/r)) 4 , where r is the microscopic time 
needed to move across one monomer. This regime appears precisely at the "stall force" 
corresponding to a disordered substrate, namely 



fs 



-In 



™7™7) /=0 



(36) 



If Ai£ has a Gaussian distribution this condition yields 

Is — r, 



(37) 



23 



The resulting behavior as the force is varied is summarized qualitatively in Fig. ^ 

It is interesting to consider the location of the stall force, f s , as well as the range of forces 
over which the displacement is anomalous, namely the region where v = limt_ >co (a;)/i = 0, 
in some more detail for both polymer translocation and molecular motors in some simple 
scenarios. These quantities characterize how the location and width of the anomalous dis- 
placement region develops as a function of temperature and chemical forces. We assume 



AE{m) with a Gaussian distribution about AE with a variance V, although it is straight 
forward to extend the results to non-Gaussian distributions with no change of the qualitative 
behavior. It is straightforward to show using Eq. (jHEJ) that the range of forces, Af, over 
which the velocity is zero satisfies 

A/ = ^ V . (38) 

For polymer translocation using Eqs. (0) and EH implies that Sinai diffusion occurs for 
the force 

/, = f , (39) 

while Eq. 0H1 implies that the range of forces around f s where the displacement is anomalous 
is given by 

Af = -L (Aj? - Ajf) . (40) 

If there are no proteins on lefthand side (cis chamber) and a small concentration, P, of 
protein is added to the righthand side (trans chamber) one can show using (|A3|) that f s ocP 
while Af oc P 2 . Thus, as the chemical bias increases both f s and Af grow. Note that in 
general one may consider proteins in both the left and right chambers. In this case even 
when the average chemical bias A/i = one may still have V > giving rise to anomalous 
dynamics even when the external bias F — 0. 

For molecular motors the situation is more interesting. The results presented above for 
the transition points between the different regimes hold even when A/i is also random. 
However, here we restrict ourselves to the simpler case when A/i is constant. In this case 
()15j) implies that for small chemical potential (A/j/T <C 1) the chemical energy difference 
AEf = o = qAfi, where q is the coefficient in the Taylor expansion of (fTKj) in A/i which is 
independent of T. Therefore, in this limit, the stall force is 

/, = ^ • (41) 
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and 



A/ 



A/i 2 



(42) 



T 



where we have assumed the purpose of a rough estimate that the chemical potential difference 
does not depend on the type of monomer. Similarly to polymer translocation, as the system is 
driven out of chemical equilibrium both f s and Af grow. However, in the limit of A/i/T 3> 1 
one obtains AEj =0 = pT where p is obtained by taking the appropriate limit in ()15|) and is 
independent of T. We then have 



implying that both quantities increase with increasing temperature. 

Note that if the force applied to the polymer or motor is held constant and the chemical 
parameters (e.g. ATP or protein concentration) are varied from their equilibrium value 
one should also observe a region of anomalous displacement (see the general expressions in 
Appendix [UJ) . These conclusions are summarized qualitatively in Figs. [TJ [2] and El 

V. EXPERIMENTAL CONSIDERATIONS 

As discussed in the previous section, the important quantity for deciding if anomalous 

2 

dynamics is present is the variance V = (AE 2 ) — (AE) of AE(m), where the overbar 
represents an average over the ensemble of random sequences. Effects related to sequence 
heterogeneity dominate when V is large compared to ksTAE. Here we estimate the ranges 
over which anomalous dynamics may be observed in experiments as well as other precondi- 
tions needed to observe this behavior. We also discuss the effect of finite time experiments on 
the shape of the velocity-force curve. In this section, we reintroduce Boltzmann's constant 



A. Polymer Translocation 

For polymer translocation, whether the variance V is large compared to ksTAE of course 
depends on a number of factors, including the base composition of the polynucleotide passing 




(43) 



and 




(44) 
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through the pore, the particular protein whose binding drives translocation, and the concen- 
tration of the binding protein. Nonetheless, it is instructive to consider an example to get 
some sense of the orders of magnitude involved. We focus on DNA binding proteins. Note 
that, like those of most such proteins, the binding sites are several nucleotides long; unlike 
in previous sections, unless stated otherwise, we will give values of V and other parameters 
normalized per nucleotide rather than per bound protein. 

The bacteriophage T4-coded gene 32 protein (gp32) is a monomeric single-stranded 
DNA ( ssDNA) binding pro t ein w hich is implicated in DNA replication and related pro- 



cesses (IColman and Oaklevl. 



1980 ). When it associa tes with ssDNA cooperatively in the 



"polynucleotide" mode ((Kowalczvkowski et all 119811 ). its net affinity 1 K nct can vary by as 
much as a factor of 10 depending on the polymer's base compositio n; in physiological sal t 
concentrations, a typical range is K ne t ~ 2 x 10 8 - 2 x 10 9 M" 1 (INewport et all Il98lh . 
In this binding mode, the binding site of each gp32 monomer is 7 nucleotides long. For 
a /xM protein concentration, K nct is large enough that almost all sites on the translocated 
ssDNA will be occupied. Upon assuming that V is determined entirely by the base depen- 
dence of K^t, we then estimate that V ~ 0.1 - 0.2(kBT) 2 for a "generic" DNA molecule in 
which each of the bases appears with roughly equal frequency. Here T is room temperature, 
ksT ~ 0.59kcal/mole. In this case, the change in free energy of a nucleotide moved from a 



buffer without any gp32 to one where the protein is present is A/i ~ k^ T (see Eq. E] 



taking the ssDNA to be a freely-jointed chain with Kuhn length 1.5 nm (jSmith et al 



Upo n 



1996). 



one finds that a force of about 10-15 pN on the polymer is required to cancel the effects of 
the protein binding. In order to have ksTAE < V, so that disorder effects can be detected, 
the value of the force must be controlled to roughly 10% or better accuracy. 



B. Molecular Motors 



To be able to measure the motion of a motor along a substrate it must remain attached 
long enough to be able to preform many moves across monomers. In other words, if the rate 
at which the motor leaves the substrate is 7 and the rate of crossing a monomer to the right 

1 The net affinity is the affinity of an additional protein molecule for a growing chain of cooperatively bound 
monomers; it differs from the affinity of an isolated protein molecule for ssDNA by an enhancement factor 
arising from the cooperative interactions 



26 



or left is or respectively, then 7 <C + w*~ must hold. In the regime of anomalous 
dynamics is of the same order of w*~. Therefore, the condition will not be fulfilled in 
this regime when the rate of hopping against the chemical bias is always very small. 
There are, however, experime nts where such a re striction does not hold. For example, 



the experiment b y Wuite et a 



( Wuite et al. 



2000) on the DNA polymerase/exonuclease 



system (see also ((Maier et all l2000l )) monitors not the displacement of a single motor but 
the location of the junction between the ssDNA and the dsDNA. Therefore, it is more 
natural to model the dynamics of the ssDNA/dsDNA junction and not of the motor. A 
motor which leaves the ssDNA/dsDNA junction is eventually replaced by a motor from the 



solution. Within ou r models this can be represen ted b y an internal state of t 



r e i un ction 



(similar in spirit to (jFisher and Kolomeiskvl . I1999T ) and (jKolomeiskv and Fishei 



mode l of this type for the DNA polymerase/exonuclease has been studied in (|Goel et al 



1999)). A 



20031 ). However, the disorder in the transition rates, present due to the heterogeneity of 



the DNA, has been neglected. In Appendix D we analyze in some detail a simple model of 
the DNA polymerase/exonuclease system. As shown in the appendix it is straight forward 
to show that the presence of heterogeneity (for example, in the energy gained from the 
hydrolysis of the different NTP's) leads to a random forcing energy landscape. One therefore 
expects a region of anomalous dynamics near where the external stretching force F' causes 
a change in direction. We stress that more realistic models with many intermediate states 
can by analyzed similarly without affecting the existence of the region with anomalous 
dynamics. Unfortunately, for this experiment an estimation of the width of the region is not 
straightforward. 

Estimates, similar to those above for polymer translocation, can be obtained for the 
random force landscapes for a number of molecular motors whi ch operate on DNA or 
RNA. Two examp 



Gelles and Landick 



e s of interest are RNA polymerases ( RNAp's) (IDayenport et al. 



1998t I.Tulicher and Bruinsmal. 



1998; 



acting on double-stranded DNA (dsDNA) (IBianco et al 



Lohman and Bjornson 



1996; 



Wang 



von Hippel and Delagouettf 



2001 



2000; 



199 8]) and he l icases 



Dohonev and Qellei 



2001 



2001). An RNAp's function is to 



transcribe DNA — that is, to synthesize an RNA "copy" with the same sequence as a DNA 
molecule. To do so, it walks along dsDNA trailing a growing RNA strand. The RNAp motor 
is powered entirely by the energy gained from the hydrolysis of successive nucleotide triphos- 
phates (NTP's) as they are added to the RNA molecule. Although the mechanism of RNAp 
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motio n is still the subject of debate ([Jiilicher and Bruinsma 



1998; 



von Hipp el and Pasman . 



20021 ) . many models suggest that at low enough NTP concentrations, its ability to move 



forward will be limited by the rate at which NTP's arrive at the catalytic site. A straight- 
forward way to force RNAp into a regime in which its motion is dominated by a random 
force energy landscape is thus to place it in a buffer with different concentrations of each of 
the four NTP's. The motor's ability to take a forward step is dependent on the incorporation 
of the appropriate NTP, and the rate of that incorporation is proportional to that NTP's 
concentration. Thus, one can in principle make V arbitrarily large and satisfy the criterion 



V > kbTAE for significant random force effects. Each factor of 10 difference between the 
concentrations of two nucleotide triphosphates, and hence in the rates to make a forward 
step, translates into a difference of kgT \n(l0) ~ 2.3k B T in AE(m). Of course, in practice 
other factors — for example the possibility that the RNAp might fall off its DNA track before 
the needed NTP arrives — will limit how large a range of concentration differences can be 
achieved experimentally. It will be interesting to see whether strong disorder effects can be 
observed. 

Another class of motors that use DNA as their track are helicases, which are needed 
to separate the two strands of dsDNA in order to facilitate various processes in the cell 
such as cell division in prokaryotes. Helicases move along the DNA by consuming energy 
from NTP's. While some helicases only b reak a few ba s e pair s at a time, others can moy e 



substantial distances along their tracks (IBianco et al 



2001 



Dohonev and Gelled . l200lh . 



Recent modelling of certain monomeric helicases (jBetterton and Jiilicher 



suggests 



that disordered DNA sequences affect helicase motion primarily through the different en- 
ergies required to open different base pairs. Random sequences thus lead to anomalous 
helic ase motion in much th e same way they do anomalous dynamics of mechanical unzip- 



ping ([Lubenskv and Nelsonl . l2000l 120021 ) . In the simplest case of "passive" opening, one finds 
that AE(m) « AE motOT + AE^^m), where AE motOT , which summarizes the forward force 
exerted by the helicase motor, is negative and has magnitude at least ~ 2k bT, and A£"dna 
is simply the thermodynamic f ree energy c o st of opening each successive base pair, with size 



roughly between 1 and 3 kT (jSantaLucia . 



1998). One thus has V ~ l{k B Tf. This large 



variance means that it should be relatively easy to observe anomalous, disorder-dominated 
dynamics in helicases as predicted earlier for DNA unzipping. If, for example, one assumes 
that the magnitude of Ai? motor is near its lower bound of 2k gT, then, in the passive opening 
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model, disorder effects should begin to appear for a mechanical load opposing the motor's 
motion of as little as 7 pN and should persist up to at least 20 pN. 

C. Finite time effects 

All calculations of quantities such as the velocity have been done by taking the limit 
of very large times and averaging over thermal realizations with the same heterogeneous 
sequence. For experiments done over finite times the velocity will not be strictly zero in 
the regime of anomalous dynamics. Instead, the velocity decays to zero as t E ~ l , where is 
the experimental averaging time used to define v as (x(£e) — x(0))/tE- The closer /i is to 
zero, the faster the decay will be. Therefore, the curve of the velocity as a function of the 
external force or chemical potential (see Figs. ^ and |2J) will be rounded becoming sharper 
and sharper as i# — > oo. To illustrate this we have carried out simulations of model (113)) 
on a single realization of the disorder averaging over thermal realizations and measured the 
v — F curve. The results are shown in Fig. El As can be seen the longer tg the closer is 
the v — F curve to that shown in Fig. ^ The convex shape of the curve near the stall force 
is clear already for averaging times tg ~ 10 5 , corresponding to motors which transverse 
distances of 0(1000) at f /T = 0. Note that, if one looks at the displacement of a single 
motor (i.e., without averaging over thermal realizations), the regime of anomalous dynamics 
will be characterized by long pauses at localized regions (corresponding to deep minima of 
the effective potential) with fast transitions between the localized regions (corresponding 
to overcoming the barrier associated with the minima). The inset of Fig. El shows a single 
trajectory of as a function of time for a given realization of disorder. The value of f /T 
was chosen to be in the region close to the anomalous velocity regime but not inside it 
(the point is at the edge of the anomalous diffusion region close to the normal diffusion 
region). As can be seen the motion of the motor is characterized by long pauses at specific 
locations along the track, with quick jumps between the pause points. The location of the 
pause points is reproducible for the same spatial disorder and different thermal realizations 
although their duration varies from simulation to simulation. Note, that since the velocity is 
finite in this regimes, over large length scales the effect of the jumps becomes unimportant. 
These pauses correspond to local minima of the effective potential and as such are inherently 
correlated with the structure of the track. Such pauses and jumps have been observed in 
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FIG. 9 The velocity as a function of f/T for different values oHe- Here A/i/T = 3 and parameters 
where chosen with equal probability to be either {p} = {5, 1, 0.3, 1, 0} or {p} = {4, 0.1, 0.7, 1, 0} (see 
text for notation). The calculated regime of anomalous velocity is 0.5116 < / < 0.699. Data was 
averaged over a 100 thermal realizations. Inset: a single trajectory shown for the same parameters 
at f/T = 0.45. 



recent experiments ((Danilowicz et all 120031 1 on DNA unzipping. 



VI. SUMMARY 



We have studied the effect of sequence heterogeneity on both polymer translocation and 
the motion of molecular motors within simple models. The models were solved exactly both 
with and without disorder. It was shown that these systems can be represented on large 
length scales and long times by a random walker moving along a random forcing energy 
landscape. Thus, in a range of forces near the stall force we expect anomalous dynamics 
where the displacement grows as a sublinear power of time. We stress again that such results 
also apply to more sophisticated models which include many internal states of the motor 
(see the discussion of the DNA polymerase / exonuclease system in Appendix D). Several 
systems in which the regime of anomalous dynamics might be wide enough to be observable 
were considered. 
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APPENDIX A: The chemical potential difference for translocating polymers 

Here we discuss the dependence of the chemical potential difference A/i for a translocating 
polymer between the righthand (trans) and lefthand (cis) sides of Fig. 0] on the protein 
concentrations and its binding energy to the polymer. Consider first a denatured polymer 
in a solution with a concentration c p of proteins which can bind to its monomers with a 
binding energy Eb < 0. We neglect cooperativity in the binding of the proteins to the 
polymer, although this effect could easily be included. Assuming an ideal solution theory 
the protein chemical potential is given \x = /Jo + Tln(P), where P = c p /c and c is the 
concentration of the solvent. Here we take the free-energy change due to an addition of one 
i solated protein to the solve nt to be /xq — T In n where n is the number of solvent molecules 



flLandau and Lifshita . Il963h . Next, we take the energy function of a polymer of length N 



inside the solution to be 

N 



H = J2(-E b cr l + fi'a l ) , (Al) 
i=i 

where cr, = 1(0) if a protein is bound (unbound) to monomer i and // is a chemical potential 
which controls the density of proteins bound to the polymer. In thermal equilibrium \i = // 
which gives for the free energy of a polymer monomer in the solution 

-Tln(l + PexpOu -£&)/T)) . (A2) 

The change in the free energy of the polymer which occurs as a result of a monomer passing 
from the left (cis) chamber to the right (trans) chamber, with ratios of protein to solvent 
concentrations Pl and Pr respectively, is given by 

Afl ~dN R ~ Tln [l + P R eM»o ~ E b )/T)J ■ (A3) 

where T is the total free energy of the polymer and Nr is the number of monomers in 
the right chamber. It is straightforward to see that this result implies that for proteins 
with different binding energy to different types of monomers A/i will depend on the type of 
monomer. 
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APPENDIX B: Deriving the effective potential from the master equation 



In this appendix we show that the equations for the probability P n (t) of being at site 
n at time t are equivalent in the long-time limit (to be specified more exactly below) to a 
random walker moving in an energy landscape constructed using E g. (HI). We demonstrate 



this 



3y eliminating the even sites from the equations of motion (see (|Lattanzi and Maritan . 
20021 ) for similar ideas). 

First, consider the equations governing the evolution of the probability, i.e. the master 
equation. For odd n one has (see Fig. EJ) 

dPnit) 



dt 

while for even n 

dP n (t) 



<Pn-l{t) + W^P n+1 (t) - « + wDP n (t) , (Bl) 



VPn-l(t) + wfP n+1 (t) - « + <)P„(i) . (B2) 



dt 

Next, we solve the equation for the odd sites and substitute into that for the even sites. 
The solution of the equation for the odd sites is 

P n {t) = e-CV-HO (p n (p) + d re (V+-D- (w^P n ^(r) + ur^iW)) , (B3) 

where P n (0) is the probability distribution at the initial time t — 0. Substituting this into 
the equation for the even sites yields 

^0 = e -(V+^)^*d Te K-+^)-(^-p n _ 2 ( T ) +wr ^P n+2 ( T )) 

+ e -K"+«a* /* rfre^+^D- + wtw^) PJr) (B4) 

JO 

- « + <)P n (t) + e-K^+^D'^p^^o) + <P n+1 (0)) . 

At times £ ^> (to^ + u>£~) one can neglect the two last terms in (jB4|) and approximate the 
integrals as follows 

' dre ( V+^K/( T ) ~ - 1 e ( V+0*/(t) (B5) 



where /(r) is assumed to vary slowly with r. In this long time approximation, Eq. (jB4|) 
reduces to 

This reduces (|B4J) to 

J t = w^w^P n+2 (t) + w^w^P n „ 2 (t) - (w^w^ + w^w^)P n (t) (B6) 

where we have rescaled times such that t — > ^j=q— ^ . As expected this equation corresponds 
to a random walker moving in a potential constructed using Eq. (j3J). 
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APPENDIX C: Derivation of the Different Dynamical Regimes 



In this appendix the expressions for the different dynamical regimes in terms of the 
hopping rates w^, w^, w^*, are given. These general equations allow a straightforward 
derivation of the expressions in the text. However, before turning to the results we outline 
the derivation of the regime where the displacement is anomalous. The derivation of the 
other regimes is much lengthier, so we only sketch the main results. 

Unless stated otherwise we assume throughout this appendix that 



log <0, (CI) 

where, as in the main text, we denote spatial averages by an overbar. Because 
exp(— AE/T) = w^wl~ /w^*w^ in our notation, this condition is equivalent to assuming 
an overall bias to the right 



AE < , (C2) 
where AE arises from the generalization of Eqs. (J^J) and (|15|) to heterogeneous systems. 



The o ther opposite regime AE > can be treated similarly. As shown by Derrida (jDerrida , 



19831 ) the velocity of a random walker on an infinite lattice model in this case is given by 



N 

lim — w , (C3) 



where 

N-l k 



En 



k=i i=i \ 



(C4) 



W i+1>i 

Here Wij is the hopping rate from site j to i. The denominator of (jC3|) can be simplified 
by replacing the sum by an average of 

1 N 

( r ) = lim T7 S Ti ■ ( C5 ) 

/V— >oo iV 7 
i=l 

Using the rates w^, w^~, w^, one finds that the average (r) is finite only if 

< 1 . (C6) 

In this case the velocity is finite. However, when the inequality is reversed (r) = oo and the 
velocity is zero. 

A much lengthier calculation along somewhat similar lines can be done to derive the other 
dynamical regimes. One obtains the following results. 
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Regime I: When 



2 



P^M <1, (C7) 
\w^w b J 

the velocity u and diffusion constant -D of the model are finite. Namely, (x) = vt and 
(x 2 ) — (x) 2 = 2Dt for long times, where the angular brackets denote an average over different 
thermal histories of the system. Assuming for simplicity that AE(m) is distributed around 



AE with a Gaussian distribution with a variance V = (AE) 2 — (AE) this condition reduces 

to 

Tl A F!\ 

> 1 , (C8) 



V 

i.e. the variance of the energy fluctuations must not be too large. Here we have used the 



fact that AE < and the relation e x = e x+ ^ x ^l 2 which holds for Gaussian distributions. 
Regime II: When 



2 



< 1 < - b - , (C9) 



the v elocity v is finite but the diffusion constant is infinite. It can be shown (jBouchaud et al 



19901 ) that in this region the long time behavior is (x) = vt and (x 2 ) — (x) 2 ~ t 2 ^ 1 , where 
1 < \i < 2. If we assume a mean value of AE with a Gaussian distribution about the mean, 

the condition reduces to 

1/2 < ' 1 < 1 . (CIO) 



V 



We have again used AE < 0. For this case it is known ([Bouchaud et all 119901 ) that the 



exponent /i is given by /i = 2T\AE\/V. 
Regime III: When 

f^]>l, (Cll) 

the velocity v is zero. More precisely (x) ~ t M where fi < 1. The diffusion about this drift 
is anomalous in the sense that (x 2 ) — (x) 2 ~ t 2>l . Assuming again a mean value of AE with 
a Gaussian distribution about the mean leads to the condition 



T\AE\ 
V 



< 1/2 , (C12) 



where again we have used the fact that AE < 0. 

Sinai diffusion: When the average bias is exactly zero, 



logl^)=0, (C13) 
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the system exhibits Sinai diffusion ((Sinail . 119821 ) with (x) = and (x 2 ) ~ (ln(t/r)) 4 , where 
t is the microscopic time needed to move one monomer. Thus, we are now considering the 
0. 



case AE 

Note that when 



log 



W a W b 
™a™b 



> 



(C14) 



namely a reversed bias where AE > 0, similar regions can be found by interchanging — > and 
<— . For example, when 

/ in 'in. \ 

< 1 , (C15) 



the velocity v and diffusion constant D of the model are finite. Such results of course require 
that the molecular motors remain attached when they reverse direction . 



Note also that the three regimes may be identified (jBouchaud et al 



1990) according to 



the parameter /x. In particular, we identify fi > 2 with regime I, 1 < fi < 2 with regime II, 
ji < 1 with regime III and \i = with Sinai diffusion. 



APPENDIX D: Simple model for the DNA polymerase/exonuclease system 

In this appendix a model of the DNA polymerase/exoneclease system is studied. It is 
shown how a more detailed microscopic model than those studied in the main text also leads 
to an effective random forcing energy landscape. However, in contrast to these models, the 
location of the transition points into the anomalous dynamics regime can not be calculated 
exactly in a straightforward manner. 

The model we consider is a simplified version of the model studied by Goel et al. 



(Goeletal 



20031 ). The model takes into account the two active sites of the motor, one 
acting as a polymerase while the other acting as an exonuclease. The system can be in one 
of five state denoted in Fig. El by (a) to (/). The figure represents only transitions which 
differ by a motion of the motor over a distance of one base. The full model along with an 
illustration of the experiment is shown in Fig. ^2 in state (a) the motor is attached to the 
ssDNA/dsDNA junction with the polymerase active site. In state (b) the motor uses the 
energy from the hydrolysis of NTP in order to be able to extend the dsDNA. States (c) and 
(d) represent similar states but now with the motor connected to the junction using the ex- 
onuclease active site. Here the motor does not utilize energy from the hydrolysis of NTP but 
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FIG. 10 The possible states of the DNA polymerase / exonuclease model. Each pair of either 
(a),(/) or (c) states differ by an addition (or removal) of one base from the dsDNA. 

instead uses the binding energy of the NMP. State (/) represents the motor unbound from 
the junction. One of the motors in the solution can bind to the junction in through either 
the polymerase or exonuclease active site. Clearly, the model is not a strictly one-dimension 
model but corresponds to a random walker moving on two lanes. 

The rates of transitions between the states are denoted in the figure. Explicit expres- 
sions similar to Eq. ED can easily be written down. The effect of the external stretching 
force F' acting on the ssDNA/dsDNA complex will cause transitions through the cycle 
(a) —?(b) -^(a) to be less favorable while transitions through the cycle (c) — -^(c) to be 
more favorable. 

To show that the energy landscape corresponding to the model in the presence of dis- 
order is indeed a random forcing energy landscap e we first calcu late the landscape for the 
homogeneous model. Using the results of Derrida (jDerridal . 119831 ) we study one cycle of the 
model (see Fig. ITUj) and calculate the ratio of the probabilities P a {ri) and P a {n + 2) of being 
in the two (a) states which differ by a translation of one base. Similarly, the effective energy 
difference between any two other sites can be calculated. With the help of Eq. (jC4|) this 
ratio can be shown to be given by 

P a (n + 2) 1 A 



P a (n) 



Wba B 



(Dl) 
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FIG. 11 The full model on the two lane lattice. The inset on the top depicts a cartoon of the 

experimental system. 




FIG. 12 The expected behavior of the velocity as a function of the external force F of the DNA 
polymerase / exonuclease system. The striped line represents the region over which the diffusion 
is expected to be anomalous. 



with 



and 



+ W^ d w7 c W }a W cf wZw2 + W^W^WfaW^WcfWfa , ( D2 ) 



B = w2w fa W cf w7 c w7 d + W^WfaWcfW^W^ + W^W^W cf W^W^ + W^W^W afW^W 



cd 
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+ W 2 W Z W 7d W 7c W cf ■ (D3) 

The effective energy landscape can be inferred by assuming an equilibrium distribution so 
that 

P (n + 2) 

p a{ ^ ' = eM(E(n) - E{n + 2))/T) , (D4) 
and the effective energy difference is given by 

AE = E(n + 2) - E(n) = —Tin ( P «( w + 2 ) ] . ( D5 ) 

V P a( n ) J 

It is now clear, using Eq. (jDl|) and arguments similar to those of Sec. IV, that if the 
set of rates becomes site-dependent, a random forcing energy landscape will develop. The 
only difference from the simple soluble models studied in the main text is that the random 
walke r representing the system is moving on a two-lane lattice. On general grounds (jFisher . 



19841 ) . this will not make a difference on the long time and large length scales behavior of 



the system. Again, one expects a region when the velocity is anomalous. The expected 
behavior of the velocity as a function of the external force F 1 is sketched in Fig. Again, 
we expect that the singularities at F^ and F^ become rounded when v is defined by a finite 
experimental time window tg, with a plateau at zero velocity becoming more and more 
pronounces as £g — > oo. 
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